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The Casimir interaction between one-dimensional metallic objects (cylinders, wires) displays unconventional 
features. Here we study the orientation dependence of this interaction by computing the Casimir energy between 
two inclined cylinders over a wide range of separations. We consider Dirichlet, Neumann and perfect metal 
boundary conditions, both at zero temperature and in the classical high temperature limit. For all types of 
boundary conditions, we find that at large distances the interaction decays slowly with distance, similarly to 
the case of parallel cylinders, and at small distances scales as the interaction of two spheres (but with different 
numerical coefficients). Our numerical results at intermediate distances agree with our analytic predictions at 
small and large separations. Experimental implications are discussed. 

PACS numbers: 45.70.-n, 45.70.Mg 



I. INTRODUCTION 

Collective phenomena in charge density and current fluctu- 
ations are known to generate Casimir-Lifshitz interactions 1 1 
2 1 that are unconventionally long-ranged for one-dimensional 
metallic objects as cylinders or wires (3JH1- These interactions 
attract currently substantial interest since Casimir forces could 
play an important role for micro and nano devices [5 1. In gen- 
eral, the building elements of such devices cannot be modeled 
as planar surfaces (as in Ref.Q]and|2]); more common shapes 
are one-dimensional structures as wires or beams. Interac- 
tions of Casimir or (non-retarded) van der Waals type between 
quasi one-dimensional shapes are also important in biological 
systems where rodlike particles auch as DNA, viruses, or mi- 
cro tuboles interact due to correlated charge fluctuations |6]|7). 
Due to correlation effects that are particularly strong for con- 
ducting objects, the interaction is not properly described by 
pairwise summation approaches. For non-conducting objects, 
these summation schemes usually yield the correct scaling of 
the interaction energy with distance and dimension of the ob- 
jects. However, for metals the scaling is different and can be 
rather sensitive to the description of the material properties 
E] . Proximity force approximations (PFA) are restricted to 
short separations ifToll but for cylindrical shapes it is impor- 
tant to study interactions also at larger separations due to their 
slow decay. 

A recently developed scattering approach can be used to 
compute the interaction between various shapes, including 
cylinders, over a large range of separations |[TTl[T2l . In the fol- 
lowing we shall employ the latter technique to study the orien- 
tation dependence of the Casimir interaction between metallic 
cylinders, both for a scalar field obeying Dirichlet or Neu- 
mann boundary conditions and for the electromagnetic field 
with perfect metal boundary conditions. The temperature T 
is assumed be either zero or in the classical limit ksT ^> he. 
An understanding of the sensitivity to orientation is important 
since the measurement of cylinder interactions can be ham- 
pered by deviations from parallelism due to external pertur- 
bations. Orientation dependence plays also a role in mixtures 
of cylindrical shapes as carbon nanotubes and of the biologi- 



cal objects mentioned before where entropic effects compete 
with the direct interaction. In the limit where the length of 
the objects is much smaller than their distance, the orienta- 
tion dependence has been studied for metallic spheroids |f]~3l . 
Here we are interested in the opposite limit where the decay 
of the interaction is expected to be much slower Bl [l4ll . The 
case of non-parallel, infinitely long cylinders has been stud- 
ied by making additivity assumptions 0. Recently, the van 
der Waals interaction between crossed cylinders, ignoring re- 
tardation effects, has been computed in the asymptotic large 
distance limit 1151 . The force between inclined cylinders has 
been derived from a dilution process for anisotropic dielec- 
tric media [16] and the results resemble pair- wise summation 
results. Advanced numerical techniques have been used to 
study the case of two perpendicular cylinders of finite length 
L (capsules) ifTTl . When L is increased, the energy is found 
to approach an L-independent value. Our numerical results 
for sufficiently small distances (where a finite L becomes less 
important) are consistent with the reported values. 

The Casimir force between crossed metal cylinders has 
been measured ifTSl . More recently, Decca et al. discussed the 
possibility to measure the thermal Casimir force and its gradi- 
ent between a plate and a microfabricated cylinder attached to 
a micromachined oscillator |[T9l . 

Before we consider the interaction at arbitrary separations, 
it is instructive to review the prediction of the PFA. This 
approximation is often applied to describe the interaction at 
surface-to-surface distances I — d — 2R much shorter than the 
radii of curvature. For parallel cylinders of length L — >• oo, 
this approximation yields at zero temperature 1(141 
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so that one has a finite energy per length in the limit of in- 
finitely long cylinders. For inclined cylinders with inclination 
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angle 9 e [0, 7r/2], cf. Fig. [I] we obtain (see App.|B|) at zero 
temperature 
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The energy is no longer extensive in L but has a simple orien- 
tation dependence with a divergence ~ 1/0 for the approach 
of the parallel configuration. In the limit of close approach, 
/ — > 0, the energy of inclined cylinders is less divergent in I 
than in the parallel setup. For inclined cylinders at short dis- 
tance one would expect a scaling of the PFA energy that re- 
sembles the one for two compact objects as, e.g., two spheres. 
Indeed, the scaling of the PFA energies of two spheres with 
I and R is the same as in Eqs. ( |I.3| l, ( |1.4) . In the following, 
we shall describe how the interaction at larger separations de- 
viates from this simple PFA estimates. We find that the in- 
teractions are no longer simple power-laws in I but acquire 
logarithmic factors and that for inclined cylinders the decay at 
very large I ^> R is even slower than for parallel cylinders. 

The rest of this work is organized as follows. In Sec. |ll] 
we describe the system and review briefly the scattering ap- 
proach that we use in the following sections and derive the so- 
called translation matrices that couple multipole moments of 
inclined cylinders. In the following Sec.|lll]we obtain analytic 
results at asymptotically large separations while in Sec. IV we 



compute numerically the interaction at intermediate distances 
and compare it to PFA predictions. Finally, we present a con- 
clusion and discussion of experimental implications in SecfV] 
Mathematical details of the derivation of translation matrices 
and the PFA formulae are provided in App. [A] and App. [B] 
respectively. 



II. SCATTERING APPROACH 

A. Geometry and Casimir energy 

Recently, a scattering approach for Casimir interactions be- 
tween non-planar objects has been developed ifTTl [121 . It 
relates the Casimir energy to the electromagnetic scattering 
properties of each object. The Casimir free energy at temper- 
ature T can be written as 



E = k B T 



OO ' 

E 

n=0 



logdet[I-N(K„)] 



(H.5) 



where the primed sum runs over Matsubara frequencies n n — 
l-KksT I [he) with the n = term weighted by 1 /2. For a sys- 
tem of two objects, this matrix is N = T1U12T2U21. Here T, 



is the T-matrix of the i object, which accounts for the geo- 
metrical and material properties of the object. are transla- 
tion matrices that describe the conversion from regular waves 




Figure 1 . (color online) Definition of distance and orientation of two 
inclined cylinders with axis-to-axis separation d and inclination an- 
gle 0. 



in the coordinate system of object i to outgoing waves in the 
system of object j. 

In the high temperature limit, fcgT ^> he, only the first 
Matsubara frequency contributes, and the energy can be writ- 
ten as 



lim E= k * r logdet [I - N(0)] = -Ml V* 1 

p=l 1 



TrpnO)], 
(H.6) 

where the relation log det A — Tr log A and the expansion of 
log(l — x) for small x have been used. The latter expression 
corresponds to a multiple scattering expansion since each fac- 
tor of N describes two scattering events, one at each of the 
two objects. The Casimir energy at zero temperature can be 
written also as a multiple scattering expansion, 
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For cylinders, the T matrices are generally known |[T2l . The 
U matrices are only known for cylindrical coordinate systems 
with parallel axes |[T2l . For non-parallel cylinders, Eqs. ( |II.6[ ) 
and ( |H.7| i remains valid. But now the U matrices depend not 
only on the distance d between the cylinder axes but also on 
their relative orientation, described by the inclination angle 
9. In particular, we will study the case of a translation along 
and rotation about the y-axis with x' = Rgx + dy where the 
matrix Rg describes a rotation about the y-axis by an angle 9, 
see Fig.[T] 
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B. Translation matrices 

In this section we derive the translation matrices U for 
scalar and electromagnetic partial waves in non-parallel cylin- 
der coordinate systems. In particular we transform outgoing 
cylindrical waves with imaginary frequency uj — in to reg- 
ular cylindrical wave in a coordinate system translated by a 
distance d and rotated by an angle 9. 



1. Scalar waves 

We consider regular and outgoing cylindrical wave func- 
tions which on the imaginary frequency axis are given by 
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with p = -J 'k 2 + k 2 and p = ^/ x 2 + y 2 . The outgoing waves 
at position x' = R^x + dy can be expanded in terms of regular 
waves at x by the linear relation 

°° T POO 
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which defines the translation matrix V n 'k' z ,nk z (d, 9). The 
mathematical details of the derivation of the translation ma- 
trix are given in App. [A] For 9 > 0, the matrix elements are 
given by 



coordinate systems. The outgoing vector waves are defined in 
terms of the outgoing scalar cylinder waves as 
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(11.14) 



and equivalently for regular vector waves. To derive the trans- 
lation matrices, we first multiply both sides of Eq. ( |II. 10| i by 
z and apply ^V'x on both sides of this equation. Then the 

left hand side gives the components of M° ut L in the coordi- 
nate frame (x' , y' , z'). To obtain the components in the frame 
(x, y, z) we multiply both sides with the inverse rotation ma- 
trix R^ 1 . This yields, using V' = RgV, 
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(11.15) 

where we have used that z = — sin(0)x + cos(#)z. This 
shows that we need to express V x (V^ E fe (x)x^j in terms of 

M^ s fc (x) and N^ g fc (x). One can show (by expressing V and 
x in cylinder coordinates) that 
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where p' = yj n 2 + k' z 2 , £ = k x /p, ? = k'Jp', k x = 
(cos(9)k z - k' z )/sm(6) and k' x = (k z - cos(6)k' z )/ sm(9). 
For 9 < the overall sign of U n >k> , n k z (d, 0) must be changed. 
The translation matrix elements for the inverse transformation 
(rotation and translation) are obtained by d — >• —d, 9 — s- —6, 
see App. [A] When we indicate that the translation matrix cor- 
responds to the inverse coordinate transformation by the argu- 
ments (— d, — 8), we have 

Un'fei.nfe, {-d, -8) = {-l) n+n 'Un>k>,nk z K 6). (11.12) 



2. Electromagnetic waves 

For the EM field one expects that the translation matrix cou- 
ples the polarizations due to the different orientation of the 



When this result is substituted into Eq. ( |II15) , we need to shift 
the index n in order to express the right hand side in terms of 
vector waves with the same index n. This can be done by 
using 



^n'K,n±i fe, (d, 9) =Ti(Z- V^Te) ' V n , K , nkz (d, 9). 

(11.17) 

This yields after some elementary algebra the translation for- 
mula for vector waves, 
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Here we have used that N = x M and x N = -M 
to obtain the translation formula for N. The translation for- 
mula for the inverse coordinate transformation is given by 
Eq. ( IIL18I) wi th Un/fc/.nk. replaced by {-\) n+n '\] n , K >nfc ., 
see Eq. ( 11.12) . Notice that the inverted sign of sin^ in the 
expression for z' is compensated by sign changes that result 
from Eqs. ^TJ7\ and pLT2] >. 
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C. Scattering amplitudes (T-matrices) 

1. Scalar field 

For Dirichlet (D) and Neumann (N) boundary conditions, 
the scattering amplitudes of a cylinder of radius R are given 
by the expressions, 
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Sn,n'S(k z - k' z ) , (11.20) 



2. Electromagnetic field 

For a perfect metal cylinder, the scattering amplitude does 
not couple electric (E) and magnetic (M) polarizations. The 
E (M) polarization is described by D (N) boundary conditions 
so that the scattering amplitude is given by 



T 



EE 
k' 



= T 



D 



T 



MM 



- T N 



(11.21) 



III. LARGE DISTANCE EXPANSION 

In this section we obtain a large distance expansion of the 
scalar and electromagnetic Casimir energies at zero and high 
temperature. The approximations performed in this section 
are valid if the cylinder radius that is small compared to the 
distance between the cylinders. To perform the expansion we 
apply the relation logdetyl = TrlogA and expand log(I — 
N) fa -N in Eqs. ( |II.6| l and $L7) . In addition to that, we 
perform an expansion of the scattering amplitudes T in the 
radius R. 



so that an effective length L e ff over which the inclined cylin- 
ders interact in the limit — s- can be determined by the 
relation sin^ = Tid/L e g. 

In order to study the high temperature limit ksT ^ hews 
need to consider only the zero Matsubara frequency contri- 
bution to the energy. A simple scaling analysis shows that 
N fe,,ofc,(K = 0) « 1/|A; Z | for k z with logarith- 
mic corrections. Hence the trace of N is not well defined, 
i.e., Nofc 2l ofc^( K = 0) is not a trace class operator, so that 
det[I — Nofc.,0 k z ( K — 0)] i s not we ^ defined. Kenneth and 
Klich showed that N is a trace class operator for compact ob- 
jects ll20l . While parallel cylinders constitute effectively a 2D 
problem of two compact discs, two tilted cylinders are non- 
compact objects whose geometry cannot be reduced to a lower 
dimensional one of compact objects. However, we expect the 
force between tilted cylinders to be well defined. This implies 
that the operator dd^ok z .vk, (k = 0) should be a trace class 
operator. This is indeed the case. In the high T limit we obtain 
the force 



Fn 
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The energy can be obtained from the force by integration, 
E T = J d °° dd' F T (d'), leading to 
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This form of Casimir interaction appears to be the one with 
the weakest decay known to date. 

For Neumann boundary conditions no logarithmic interac- 
tion appears but the energy is proportional to the product of 
the cross-sectional areas of the cylinders. In Eqs. ( |II.6[ ) and 
( |H.7| i we have to consider the terms with p = 1, \n\ < 1 and 
|n'| < 1 since I' n (z)/K' n (z) - + O [z 4 ] for \n\ < 1. 
With this expansion we get to lowest order in R/d the energy 
at zero temperature, 



A. Scalar field 

To obtain the asymptotic Casimir energy for Dirichlet 
boundary conditions at large distances d ^S> R, we need to 
consider only the terms with n = nf = and p = 1 
of Eqs. ( |H.6| l and ( |H.7| i. Using for the small radius expan- 
sion of the scattering amplitude T that for small z one has 
Io(z)/K (z) w — log _1 (z), we get for the energy at zero 
temperature 
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The corresponding energy for parallel cylinders of length 
L -> oo is ifTHl 
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which can be compared to the corresponding energy for par- 
allel cylinders |[T4ll . 
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so that the effective length L e s in the limit 8 — » is deter- 
mined by sin 9 = (3ir/8)d/L e g. In the high temperature limit 
the energy can be obtained directly since N(k — 0) is a trace 
class operator, yielding 
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In contrast to the Dirichlet case, the energy has a more com- 
plicated orientation dependence. 
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B. Electromagnetic field 



We calculate the energies with the help of the expansions 
given in Eqs. (|II.6|l, (|H.7|i. In general, the translation ma- 



trix of Eq. ( |II. 18 1 mixes the two polarizations. However, the 
leading part of the energy at far distance is given by the E- 
polarized waves only due the logarithmic dependence on fre- 
quency of the scattering amplitude for E modes with n = 
and a power-law dependence for M modes. In the high tem- 



perature limit where all matrix elements are computed in the 
limit k — y 0, the polarizations are decoupled at all distances 
since the matrix elements that couple different polarizations 
are proportional to k, see Eq. ( |II.18| l. We find that the asymp- 
totic Casimir energy for perfect metal cylinders at zero tem- 
perature can be written as 
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where the function 0(9) is defined as 
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The orientation dependence described by 0(9) results from 
the geometric factor cos 9 — (k z £/p) sin 9 that appears in the 
coupling of the E polarization in Eq. ( |II.18[ ). We have O(0) = 
1 so that for 9 — > the Dirichlet result given in Eq. ( |HI.22| > 
is recovered. For perpendicular cylinders one has 0(ir/2) = 
1 — log(2). At intermediate angles the function 0,(9) can be 
computed by numerical integration. The corresponding result 
is shown in Fig. [2] The function 0(9) can be also expanded 
as a Fourier series 0(9) = X^^Lo cos (2n#) with the low 
order coefficients given by Slo = 0.6137, O2 — 0.3333, O4 = 
0.0333 and 6 = 0.0096. 
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Figure 2. Amplitude function fi(#) of the asymptotic electromag- 
netic Casimir energy of Eq. l |III.30| l as function of inclination angle 



Next, we consider the high temperature limit. In this limit, 
at all separations the polarizations are decoupled since the 
coupling is proportional to k. This implies that the interaction 
is the sum of the energy for Dirichlet and Neumann bound- 
ary conditions. However, the presence of a (n-independent) 
geometric factor in the translation matrix elements for both E 
and M polarizations, see Eq. ( |H.18| l, could modify the inter- 
action in the electromagnetic case. That this is not the case 



r 



can be seen by considering the limit k — >• of the geometric 
factor (p/p')[cos9 — (k z ^/p) sin9] which is unity. Hence, in 
the high temperature limit, the interaction is precisely the sum 
of the Dirichlet and Neumann energy. At large separations, 
the Dirichlet contribution dominates. This leads for the elec- 
tromagnetic case to the same problem as for a Dirichlet scalar 
field due to the fact that No fc z ,o k„ at k = is not a trace class 
operator. 

However, as in the scalar case, we are able to obtain the high 
temperature limit of the force, because c^No^.ofc^ (n = 0) is 
a trace class operator. The far distance approximation of the 
force is 
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4dsin(6») log z (R/d) 
and the corresponding energy is obtained by integration as 
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(11132) 



which is identical to the Dirichlet result of Eq. ( |HI.25| l. It 
should be mentioned that this slow decay with distance is a 
consequence of the metallic response of the cylinders. For 
non-conducting cylinders a power-law decay is expected. 



IV. NUMERICAL RESULTS FOR SMALLER 
SEPARATIONS 

In this Section we present numerical results for the interac- 
tion between inclined cylinders. This allows us to go beyond 
the limit of asymptotically large separations. We focus on the 
zero temperature with the energy given by Eq. (II.7) . The 
broken translational symmetry of inclined cylinders leads to 
a matrix N that is non-diagonal in the continuous index k z . In 
contrast to compact objects or systems with translational sym- 
metry (parallel cylinders), one has to compute the determinant 
of a non-diagonal operator over a continuous index. There 
are basically two approaches to obtain an approximation for 
the determinant. One can approximate the determinant by the 



6 



trace of a series of powers of N [see Eq. jll.7\ ] which is known 
as multiple scattering expansion or perform a discretization in 
k z and compute the determinant over the resulting finite set of 
discrete variables. While the first approach is limited to suf- 
ficiently large separations since the multiple scattering series 
needs to be truncated at some order, the second approach re- 
quires enough discretization points in order to minimize the 
discretization error. We have decided to follow the second 
route. This discussion shows that the evaluation of the in- 
teraction between inclined cylinders is substantially more ex- 
pensive in terms of numerical operations than that of paral- 
lel cylinders. The second discrete index n of the matrix N is 
truncated at some value n max as in the case of compact objects 
ifTTl . The choice of n max depends on the range of distances for 
which one would like to compute the energy. n max increases 
with decreasing minimal distance. For our results shown be- 
low, we found it sufficient to consider |n| < 3 = n max - We 
studied two different values for the inclination angle, 8 = it/2 
and 9 = tt/4, for separations d > 2.22R. 

Our numerical results for the electromagnetic Casimir ener- 
gies are shown in Figs. [3] and [4] (dots and solid curves). They 
are scaled by PFA energy for cylinders with 9 = ir/2 and 
with 9 = tt/4:, respectively. Regarding the approach of the 
numerical results to the asymptotic expression for the energy 
of Eq. ( |III.29| l (dashed curves), we observe similar behavior 
as in the case of parallel cylinders lfl4ll : Since the energy 
decays logarithmically (with sub-leading logarithmic correc- 
tions), the actual energy converges to the asymptotic result 
only at extremely large separations that are not shown in the 
plots. At intermediate distances, higher order partial waves 
need to be included and are expected to lead to an approach 
of the energy to the PFA result at small separations. This 
tendency indeed can be observed for our numerical data. It 
should be noted that beyond the distance d « 3.33i? where 
the curve for the numerical result starts to bend downwards, 
the result becomes inaccurate and more partial waves must be 
included. Instead of performing computations with more par- 
tial waves, we compare our numerical results to the prediction 
of a recently developed gradient expansion of the interaction 
energy for gently curved surfaces ll2Tll22i . The latter approach 
yields for the first correction to PFA energy the result 
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where I = d — 2R is the surface-to-surface distance. This 
result is shown as dotted curves in Figs. [3] and [4] It can be ob- 
served that our numerical results nicely approach the predic- 
tion of Eq. ( |IV.33| > for short distances. Hence, our numerical 
results together with Eq. ( |IV.33[ > provide the overall behavior 
of the Casimir interaction between inclined cylinders. 

We are also interested in the orientation dependence of the 
energy at a fixed distance. As can be seen from Eq. ( |III.29| > 
and Eq. ( |I.3[ ), the energies have a different angular depen- 
dence. To study the angular dependence of the energy at in- 
termediate distances, we show in Fig. [5] the rescaled energy 
ui(r, 9) — Eo(r, 8) sin(6*) as function of the inclination angle 
compared to the same function at 8 = it/ 2, for different dis- 
tances. One observes that at large distances the dependence 
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Figure 3. (color online) Electromagnetic Casimir energy normalized 
to the PFA energy as function of the inverse distance for perpendic- 
ular cylinders (9 = n/2). The dashed curve represents the asymp- 
totic energy of Eq. l |III.29[ l and the dotted curve shows the energy of 
Eq. |TV33]|. 
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Figure 4. (color online) Same plots as in Fig. [3] but for inclination 

9 = tt/4 



of cj(r, 8) on 9 is close to the asymptotic result of Eq. (111.29 1, 
while for reduced distances w(r, 9) becomes more flat, indi- 
cating that it converges to the constant value of the PFA result. 



V. SUMMARY AND CONCLUSION 

We have studied the interaction of inclined cylinders with 
Dirichlet, Neumann and perfect conductor boundary condi- 
tions. In order to study this system, we have obtained the 
translation matrices between outgoing and regular cylindrical 
waves in inclined coordinate systems, both for scalar and elec- 
tromagnetic waves. The interaction energies have been com- 
puted over wide range of separations, using analytic and nu- 
meric approaches. The zero temperature and classical high 
temperature limits were considered. In the latter, we have 
found that the zeroth order Matsubara term corresponds to 
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Figure 5. (color online) Ratio cj(6) /u>(n /2) with 10(6) = 
Eo(r, 6) sin(#) and r — R/d as function of the inclination angle 
8. The dashed curve represents the asymptotic result of Eq. l |lil,29} . 
The solid curves connect numerical data for r = 0.01, r — 0.1, 
r = 0.2, and r — 0.3 (from top to bottom). The PFA yields for the 
ratio unity. 



a non-trace class operator N(n — 0) for a scalar field with 
Dirichlet boundary conditions and for the electromagnetic 
field. However, we could obtain the energy by integration 
of the force which is well-defined since ddN(K — 0) is a 
trace-class operator. We should remark that the zeroth order 
Matsubara term contributes to Casimir energy at any non zero 
temperature so that this result is relevant not only in the high 
temperature limit. The non-trace class property of N(n = 0) 
can be related to the metallic response of the cylinders. As 
soon as one of the two cylinders is non-conducting or one of 
the cylinders has a scalar field boundary condition different 
from the Dirichlet case, N(k = 0) is a trace class operator. 

We note that our analysis for perfect metal cylinders could 
be easily extended to dielectric cylinders since the translation 
matrices remain unchanged. In this case it is expected that at 
large distances the energy decays according to a power-law, 
the zero Matsubara frequency contribution leads to a trace 
class operator, and the interaction is generally weaker than 
in the perfect metal case. 

The Casimir interactions between cylinders as a function 
of their inclination angle 9 could be experimentally probed. 
The force between cylinders with fixed inclination could be 
measured directly or by measuring the effect of the Casimir 
force on the mechanical response (oscillations) of thin metal- 
lic wires. The rather slow decay of the interaction between 
inclined cylinders might allow to consider separations that are 
well beyond the validity range of the PFA. Another experi- 
mentally interesting quantity is the torque t = —dE/dQ be- 
tween the cylinders for which a torsion pendulum might be 
employed. For future work, it would be interesting to study 
the influence of finite conductivity and the related anomalous 
scaling of the energy [9] on the orientation dependence. 
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Appendix A: Derivation of translation matrices 

In this appendix, we provide the derivation of the result in 
Eq. (JTLTTjl. We start from the Fourier transform of outgoing 
waves in the xy-plane, 



K n {pp)t 



k-r. " 



OO u 
-OO 27T 



2tt 



p 



v » 6 (A.l) 



k 2 +k 2 +p 2 ' 



The integration over k y can be easily performed using the 
residue theorem. Since the integrand has poles at k y — 

±iy/p 2 + k 2 , we obtain 



K n {pp)e me = {-%)" 



k x T \Ax + P 2 



Jk^x^y^Jkl+p 2 

2^kJTp^ ' 



(A.2) 



where the minus (plus) sign applies to y > (y < 0). With 
the coordinate transformation x' = Rgx + dy and the corre- 
sponding wave vector components k x — k' x cos 9 — k' z sin 9, 
ky = k y = ±i\J k 2 + k x + k 2 , k z — k' x sin 9 + k' z cos 9 we 
get the outgoing wave in the primed coordinate system for 

y' > o (V < 0) 



/ k' x =F \ k' x 

dk' 

OO \ P 



C? fei 0O = H)' 

e i(k m x+k y y+k z z) e T^/K?+p^ 
X 



2 + P ' 2 



2xlk' 2 +p' 2 



(A.3) 



with p' = y k 2 + k' 2 . Now we use the expansion of 2D plane 
waves in cylindrical waves, 



e i(xk z +k y y) = ^ % n j n (p^jkl + kfj e me 



-in arccos 



xe 



(A.4) 



With the relations ip = \jk^-\- and 



-in arccos 



k x -\- iky 

p 



(A.5) 
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the plane wave can be written as 



nGZ 



(A.6) 



By insertion of this sum in Eq. ( |A. 3| >, we can express the out- 
going wave in terms of incoming waves, 



(pp) e in V fc ** , (A.7) 



where we assumed y 1 > 0, rf > 0, and > 0. In order to bring the result into the form of Eq. (LLIO) , we change the variable of 
integration from k' x to k z using k z — cos(9)k' z + sin(0)fc^, and dk' x — dk z / sin 6*, which yields 



/oo 
dk 



/ / \ " - 

(_ z )n+n' lk' x -^k'l+ p' 2 \ U _ ^fc2+p2\ " e -<V*'2+P' 2 



sin(0) 



2Jk>'+p" 



r 



where we have used c/>™f (x) 
understood that 



= L 



(pp) e m9 e tkzZ , and it is 



k x — k z cot 9 



sin t 



k z 



sine 



- ^ cot . (A.9) 
fe^/V Eq. dAj) is 



With the definitions £ = fc^/p and £' 
equivalent to Eq. ( |II.l 

For 9 < the overall sign of U n /fc< )n fe a (ci, 0) must be 
changed. The translation matrix elements for the inverse of 
the transformation x' = Rgx + dy is obtained by assuming 
y' < and d — > —d, 9 — > —9. When we indicate by the argu- 
ments (— d, —9) that the translation matrix corresponds to the 
inverse coordinate transformation, we have 



U, 



(-l) n+n U n , KM ,(d,6). (A.10) 



Finally, we note that the translation matrices for parallel cylin- 
drical coordinate systems follow from a direct integration of 
Eq. dA7) over k' x in the limit 9 — > 0. 



Appendix B: Proximity Force Approximation 

In this appendix we derive the PFA for the energy of two in- 
clined cylinders at zero temperature, Eq. ( |I.3| l, and in the high 
temperature limit, Eq. ( |I.4[ ). The area across which the two 
inclined cylinders overlap, viewed along the axis that is per- 
pendicular to the two cylinder axes and that intersects the axes 
in their crossing point, forms a parallelogram of edge length 
2i?/sin(9. Let us denote the coordinates along the edges of 
this parallelogram as u and v. Then the local distance h(u, v) 
between the two cylinder surfaces, measured normal to the 
plane that is spanned by the cylinder axes, is given by the 
function 

h(u, v)=d- \Jn 2 - (usinfl-i?) 2 



R 2 - (vsm9 - R)' 



(B.l) 



where d is the distance between the cylinder axes. Taking into 
account that a surface element of the parallelogram is given 
by sin(0) dudv, the PFA energy at zero temperature can be 
written as 



PFA 
0,6 



hcTT 2 

720 



2R 



ffi^ sin(9)dudv 
h 3 (u, v) 



(B.2) 



When we introduce the surface-to-surface distance I = d — 
2R, after performing a change of integration variables the en- 
ergy can be written as 



pPFA 
^0,6 



s/R/l 



IT 2 he 1 I 
720 sin0 R? 

y/R/l 



dsdt 



fR/l 



i-J^-V 1 -^ 2 



(B.3) 



This expression has the advantage that we can expand the 
square roots for small I / R which leads to 



E> 



PFA 



ir 2 he 1 R 
720 sin(0)Z2 



dsdt 



.(B.4) 



Wl J-yJW [1 + 5 {s 2 + t 2 )} ' 

In this expression we can extend the integration limits to infin- 
ity to obtain the limiting behavior for small I / R. This yields 



,,_ pPFA _ T he 1 R 

<™ <* " 720 M0) I 2 Jo 



ir 3 he 1 R 



pdp 



720 sin(0) I 2 



(B.5) 
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which is the result of Eq. ( |I.3[ ). For small I / R the latter ap- 
proximation deviates from the exact integral of Eq. ( |B.3| l by 
less than 1% for l/R < 0.01. It is instructive to compare 
this PFA energy to the one for two spheres of radius R and 
surface-to-surface distance I which is 



f 11 



PFA 
spheres 



7T 3 HcR 

1440 P 



(B.6) 



and hence reduced by a factor of 1/2 compared to the case of 
perpendicular cylinders with 8 = tt/2. 

In analogy to the above computation, it is possible to ob- 
tain the PFA energy for the high temperature limit. Using the 
high temperature interaction between two perfectly reflecting 
plates, the PFA energy for two inclined cylinders at high tem- 



peratures is given by 

C(3) 



&CL9 



k R T 



8tt 



in 



in 

■h,:fH 



sin(0)dudv 
h 2 (u, v) 



(B.7) 



Performing the same transformations and approximation as 
for T = 0, we find for small / /R the PFA energy 



lim E 



PFA 



= -k B T 



-kpiT 



C(3) 1 
8?r sin(0) 

C(3) 



1 R 



2tt 



dip 



pdp 



4 sin(0) 
which corresponds to Eq. \\A\. 



(B.8) 
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